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mal distributed. Multi-dimensional central limit theorem is also ob¬ 
tained here. These results have several applications to various physical 
models and random matrix models, such as the Anderson model, the 
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1 Introduction 


We are here concerned with the tridiagonal random matrix 


/ di bi 


\ 


Qn 


Oi (^2 

0-2 ds bs 


( 1 . 1 ) 



where a*, bi and di, 1 < i < n, are random variables with oq = &n = 0. 

Tridiagonal random matrix attracts signihcant interests in various helds. 
In quantum mechanics, it is a hnite-difference model of one-dimensional ran¬ 
dom Schrodinger operator, such as the extensively studied Anderson model 
where a* = 6^ = — 1 and all di are independent and identically distributed 
(see e.g. [281 [7]). The non-symmetric model with tti/bi > 0 also arises in 
the non-Hermitian quantum mechanics of Hatano and Nelson, see e.g. [21] 
and references therein. Tridiagonal random matrix is also a basic model of 
random walks with random environment in chain graphs, interesting exam¬ 
ples include the random birth-death Markov kernel proposed in [5] [8] , where 
tti-i + di + bi = 1 and {(aj_i, di, b^)} is an ergodic random held, and the ran¬ 
dom birth-death Q matrix recently studied in [23], where di = — (ai_i -|- bi) 
and {oj} and {6^} are two sequences of strictly stationary ergodic processes. 

Moreover, tridiagonal random matrix also plays an important role in ran¬ 
dom matrix theory. One well-known model is the /3-Hermite ensemble, in 
which tti = bi, tti is distributed as iXip is the y distribution with 

i/5 degrees of freedom), (3 > 0, di is normally distributed as N{0,2/(3), and 
{oj, di} are independent. This model was hrst proposed in [m and gener¬ 
alized the classical Gaussian ensembles for f3 = 1,2,4, corresponding to the 
Gaussian orthogonal ensemble (GOE), Gaussian unitary ensemble (GUE) 
and Gaussian symplectic ensemble (GSE) respectively. Due to its simple 
tridiagonal structure, one expects to investigate some interesting spectral 
phenomena of general random matrices from the study of the tridiagonal 
random matrix fll.ll) . For the generalization of the /5-Hermite ensemble to 
general symmetric tridiagonal random matrix (with independent entries), we 
refer the interested reader to ra. which is another motivation for the present 
work. 
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The interest of this paper lies in the Gaussian fluctuations and deviations 
of the trace of powers, or in other words, the moments of empirical spectral 
distribution of the tridiagonal random matrix fll.ip . 

The fluctuations of the traces of random matrices are extensively stud¬ 
ied in the literature and, in a general context, turn out to be Gaussian. 
For instance, the traces of the classical (unitary, orthogonal and symplectic) 
compact groups were proved by Diaconis and Shahshahani jT3], by using 
the representation theory, to be independent and normally distributed in the 
limit. For the Wigner matrices, the Gaussian fluctuations of the traces are 
presented in [H Theorem 2.1.31] (see also [22l Theorem 2.1]), based on elab¬ 
orate computations of moments. See also [22] for the symmetric tridiagonal 
random matrix with indepenent entries, where the Gaussian fluctuations of 
the traces were obtained by a judicious counting of levels of paths. We also 
refer to [21 [181 iMl iSSl [301 El] and references therein for fluctuations of other 
linear eigenvalue statistics. 

Furthermore, large deviation was well-known for the empirical spectral 
distribution of the unitary invariant ensembles, including the /9-Hermite en¬ 
semble (see [3] and DP)- See also [1] for the non self-ajoint matrix with 
independent Gaussian entries. Moderate deviation principle for the empiri¬ 
cal spectral distribution of the Gaussian divisible matrix was studied in [T2] . 
Moreover, for the eigenvalue counting function of the determinantal point 
process and Wigner matrix, we refer to [U] and [TB] . 

Motivated by the physical models and random matrix models mentioned 
above, we here consider the tridiagonal random matrix fll.ip under quite gen¬ 
eral assumptions (see {H.l) and {H.2) in Section [3l see also (Ff.3) in Section 
up, which, in particular, allow the non-independent entries in fll.ip . We prove 
the Gaussian fluctuations of the traces of such model. Multi-dimensional cen¬ 
tral limit theorem is also given. Moreover, the large deviations of the traces, 
in relation to those of additive functionals of uniform Markov chains, as 
well as the moderate deviation results are also obtained here. These results 
are applicable to various models, such as the Anderson model, the random 
birth-death Markov kernel, the random birth-death Q matrix as well as the 
symmetric tridiagonal random matrix, including the /3-Hermite ensemble. 

Our proof is quite different from the standard method of moments men¬ 
tioned above and relies crucially on a path expansion of the trace (see fl2.5p 
below), which is formulated according to the types of circuits determined by 
the tridiagonal structure of fll.ip . This formula was recently used in our arti- 
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cle 123] to study the limiting spectral distribution of the random birth-death 
Q matrix. The advantage of this formula is that, it reveals that the fluctu¬ 
ations and deviations of the traces are asymptotically the same to those of 
the sum of m-dependent random variables. This new point of view gives a 
unihed way to understand the limit behavior of traces of quite general tridi¬ 
agonal random matrices with non-symmetric structure or non-independent 
entries. Moreover, it also allows to employ the analytic tools (e.g. blocking 
arguments, Lyapunov’s central limit theorem and large deviations of addi¬ 
tive functionals of Markov chains) to obtain the Gaussian fluctuations and 
deviations of the traces. 

The remainder of this paper is organized as follows. In Section [2] we hrst 
set up some preliminary notations and dehnitions, then we present the path 
expansion formula fl2.5p . Section |3] is devoted to the Gaussian fluctuations 
of the traces, and Section 0] is concerned with the deviation results. Finally, 
Section |5] includes some discussions on the main results of this paper, and 
the Appendix, i.e. Section [3] contains some technical proofs. 

Notations. Throughout this paper, for x G M+, [x] denotes the largest 
integer not greater than x, f = 0(g) means that \ f /g\ stays bounded, and 
fn = o(l) means that |/„| tends to zero, as n —)■ oo. 

2 Preliminaries 

Let us hrst recall some notations from our recent paper [23]. In particular, 
we will classify the circuits according to their types. Then, we shall formulate 
the path expansion formula of the trace precisely. 

Let k >1 he hxed. For every 0 < / < [|], 1 < i < u — set 

i-i i 

Qw"’*' (2-1) 

j=0 j=0 

where, ut; = (mi, m 2 ,..., m^) and iti = (uq, ui,..., u;). We say that r^i and 
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( 2 . 2 ) 


Ui are admissible, if rrij > 0 and > 0 for every 0<h<l, 

i i 

2 ^ rrij + ^ rih = k, 

j=l h=0 


and if for some I > 0, there exists 1 < p < /, such that = 0, then rrij = 0 
and = 0 for all p < j < /, p < h < 1. Set 


\T;. 




The intuitions of these quantities can be seen as follows. By the tridiag¬ 
onal structure of fll.ll) . we note that 


k 

TrQi = Y, ( 2 - 4 ) 

TrS^n t = l 

where denotes the set of all circuits, i.e. 


K = {tt : (I,-- - ,fc) ^ (I,-- - ,77.) : \Tij -TTj+il < 1,1 < j < A;,7rfc+i = tti}. 

The types of can be determined by the set and the vertices i, 1 < i < n. 
Indeed, I is the largest length between various vertices in the circuit tt, and i 
is the leftmost vertex in tt. The vector r^i determines the circuit n consisting 
of the subindices of off-diagonal entries in Q.,^, and the coordinate mj+i is 
half of the number of edges with the vertices i + j and i + j + 1. Similarly, 
the vector iti is related to the diagonal entries in and the coordinate nj 
is the number of loops with the vertex i + j. Thus, dehned in fl2.ip 

can be viewed as a representative element of the type and the 

leftmost point i. 

Let denote the number of circuits of the same type 

and the leftmost vertex i. Note that for every 1 < i,j < n, Ci^' * = 

Cj^ji’^i. That is, the number of circuits with the same type {/, ht;, itj}, 
though different leftmost vertices, are also the same. Hence, we can set 

With these notations, the expansion of the trace (12.41) can be reformulated 
according to the types of the circuits i-e. 

n—l 

TrQi= Y. Cf"^‘ ( 2 . 5 ) 
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This formula was recently used in our article [23] for the study of limiting 
spectral distribution of the random birth-death Q matrix and, as we shall 
see later, is crucial for the formulations and proofs of the main results in this 
paper. 

We conclude this section by taking TrQ\ and TrQ‘^ for examples. For 
TrQ^, in this case / = 0,1, contains the types (0, (0), (3)), (1, (1), (1,0)) 
and (1, (1), (0,1)), and is equal to 1, 3, 3 respectively. Hence, by fl2.5p 

we have 

n n—1 n—1 

TrQl = + 3 '^{aibi)di 3 '^{aibi)di+i. 

2=1 2=1 2=1 

For TrQ^, / = 0,1, 2, have the types (0, (0), (4)), (1, (1), (2, 0)), (1, (1), (1,1)), 
(1, (1), (0, 2)), (1, (2), (0, 0)) and (1, (1,1), (0, 0, 0)), and the corresponding 
1^4^ 4^ 4^ 2 and 4 respectively. Thus, it follows that 

n n—1 n—2 

aibiiAd?^ -\- Adidi^i + Ad'^j^^ + ‘2aibi) + A aj6jaj+i6j+i. 

2 = 1 2=1 2 = 1 

3 Gaussian fluctuations 

In this section, the main results Theorem 13.21 and 13.31 are formulated in Sub¬ 
section 13.11 and then they are proved in Subsection 13.21 and 13.31 respectively. 
Moreover, several applications are also given in Subsection 13.41 

3.1 Main results 

Let m be a fixed nonnegative integer. The random variables {W} are said 
to be m-dependent, if for any positive integers i and j with j — i > m, Xj 
is independent of the a-held generated by {Xh, I < h < i}. In particular, 
0-dependence is equivalent to independence. 

Motivated by the physical models and random matrix models mentioned 
in Section [H we introduce the assumptions {H.l) and {H.2) below. 

(H.l) In the symmetric case (i.e. a* = bi for all i > 1), the off-diagonal 
entries {oj} are independent, and the diagonal entries {di} satisfy that 
di = /(aj_i,aj) with / a continuous function on or {ai,di} are all 
independent. 
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In the non-symmetric case, the random vectors {(aj_i, di, 6j)} are inde¬ 
pendent. 

As regards the asymptotic behavior of the entries in fll.ip . we assume that 

(H.2) Let a and e be two nonnegative constants, 0 < e < a. There exist 
constants a, d, b and random variables rji, Q, ^i, rj, (, such that 

r"(ai_i, di, bi) = (a, d, b) + 0, Ci), (3-1) 

where r/j-i, 0 and satisfy 

^ (3-2) 

“ ” means convergence in distribution, and for each /c > 1, 


snpE(r;f+ Cf+ ef)<oo. (3.3) 

i>l 

Remark 3.1 (*). By Holder’s inequality, fl3.3p implies that all moments of 
VnXn,fn are finite. 

{ii). Assumption {H.l) allows to treat the tridiagonal random matrix 
(EB with non-independent entries. The conditions on di in the symmetric 
case is mainly motivated by the random birth-death Markov kernel and the 
random birth-death Q matrix, where /(aj_i, Oj) is of the form 1 — (ai_i -|- Oj) 
and — (aj_i -|- Oj) respectively. 

{Hi). In the case that the entries of (ED are independent and identically 
distributed (i.i.d.) and have all moments finite, the assumptions {H.l) and 
{H.2) are obviously verified. 

{iv). In Assumption {H.2) with a > 0, when 0 < e < a (resp. e = 0), 
the entries convergence in distribution to degenerate (resp. non-degenerate) 
random variables. The degenerate case actually includes the fi-Hermite en¬ 
semble. See Subsection \3.4\ for more details. 

The main results in this section are formulated below. Taking into ac¬ 
count the complicated formulations of covariances in the case £ > 0, we shall 
consider the case £ = 0 and 0 < e < a in Theorem 13.21 and 13.31 respectively, 
to make the structure of the variance and covariances clean. 
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Theorem 3.2 Assume {H.l) and {H.2) with £ = 0. 

(i). For each k > 1, let ruk = [|] (resp. [|] + 1 j in the non-symmetric 

(resp. symmetric) case. Set TrQ^ := TrQ^ — KTrQ^. Then, 

4 jv(o, Dk). (3.4) 


Here 


Dk = 


mk 


2ak + 1 
and for each 1 < i < mk + 1, 


Var{Zk,i) + 2 ^ Cov{Zk,i, Zk,i+j) 

j=i 


(3.5) 


i-i 


= E ‘ n( 5 ‘+i 5 


j=0 


(3.6) 


where pi, Q, fi have the same distributions as those of a + p,d + ( and b + ^ 
respectively, and they satisfy Assumption {H.l) with all ai,di,bi replaced by 
rjiXi and fi respectively. 

{ii) Given ki, ■ ■ ■ ,kr > l,r > 1, set mij = max{mfe .,We have 

n-h (n-'^^^TrQ’)f, • • • , ( 3 . 7 ) 


where 4>a is a r-dimensional normal distribution with mean zero and covari¬ 
ance matrix A, for 1 < j < r, 


A(i, j) 


Cov{Zk^^i, Zkj^i 


(y(^hi kj^ “ 1 “ 1 

mij 

+ {Cov{Zki^i, Zkj^i+h) + Cov{Zki^i+h, Zk,jX) 


h=l 


, (3.8) 


where Zki^h and Zk^^h, 1 < ^ + 1, are defined as in fl3.6p . 


Theorem 3.3 Assume {H.l) and {H.2) with 0 < £ < a. Let TrQ^,mk and 
mij be as in Theorem \3.2[ 

(i). For each k > 1, 0 < j < mk, there exists the limit 


crfc(l+j, 1) := limn ^'>Cov{Xk,n, Xk,n-j), 

n^oo 


(3.9) 












and we have 


n-(“"+^-^)TrQ^4iV(0,A), 


where 


Du — 


1 


ruk 


;(c’‘fc(l51) + 2 (Tfc(l + j, 1)). 
i=i 


2ak + 1 — 2e 

{ii). Given fci, • • • ,kr > 1, r > 1. For every 1 < i, j < r and 0 < h < 
niij, there exists the limit 

+ h, 1) := lim (3.10) 

n—^oo 

and we have 

n-"^+%n-^^^TrQi\ • • • 4 $a, 

where $a is the r-dimensional normal distribution as in Theorem \S.2. but 
with the covariances A{i,j) defined by 

1 


A{iJ) 


^ki,k^ (1) 1) 


(y(^ki “h kj^ 1 — 2^ 

rriij 

+ 1 + /i) + + h, 1)) 


h=l 


, (3.11) 


where (Tki,kfil, 1 + h) = akj,kiil + h,l). 

Remark 3.4 Actually, in the case 0 < e < a, ak{l+j, 1) and (l + /i, 1) 
can be calculated explicitly in terms of the covariances ofrjnyfn o-nd Cn- Since 
the formulations are complicated, we omit them in the statement of Theorem 
13.31 Concrete calculations are shown in Corollarv \3.& below for the symmetric 
tridiagonal random matrix motivated by f29f and the fi-Hermite ensemble. 


3.2 Proof of Theorem 13.2 


The key observation for the proof is that, the path expansion formula of the 
trace fl2.5p indicates that the fluctuation of TrQ^ is approximately the same 
as that of 




C. 




n 

i=l 


n 



2 = 1 
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where 




k^i 




l Q^u'^I 


{l,r^l,lti)£^h 


c, 




1-1 


n(o.+j6.+;)”'« 

j=0 j=0 


(3.12) 


{^k,i}i>i is mfc-dependent with rrik dehned as in Theorem 13.21 due to As¬ 
sumption {H.l) and the hnite width of band in the tridiagonal random matrix 
(HU. With this new point of view, we can expect the Gaussian fluctuation 
of TrQ^, inspired by the fact that m-dependent random variables are ap¬ 
proximated normally distributed in the stationary case (see e.g. [Ml El 19]). 
On the technical level, in order to deal with the non-stationary case a > 0 
in Assumption {H.2), we will employ the standard blocking arguments to 
separate the sum 'Yl'i=i ^k,i into independent blocks and small ones. 


Proof of Theorem \3.^ {i). We hrst note that 


E 


n 


-(ofc+i 


I'UTvQi-Y.Xt.i 


i=l 


0 , 


(3.13) 


where j — j. (See the Appendix for the proof.) 

Thus, we only need to consider the fluctuation of Yl'i=i which is 
actually the sum of mfc-dependent random variables. 

Let K < |, n' = p = [^], r = n — pn'. Set For 1 < 

^ in'—mis ^ ^ in' _ 

i < p, let Un,i := XI Xj- For 1 < i < p - 1, set = X Xj, 

j={i—l)n'+l j=m'—m^+l 

^ n ^ ^ p ^ ^ 

^n,p = X X,i and X = X Zn,i- Note that, {Un,i}i<i<p are indepen- 

j=pn'—mk-\-l "^=1 

dent, and so are due to the m^-dependence of {X^, j}. Moreover, 


n 


-{ak+P 


IV IV t' 

Xk,i = n~^ Yn,i = n~^ Un,i + n~^fn 


2=1 


2=1 


2=1 


Let us show that, as n ^ cxo. 


n 2 T. 


0 . 


(3.14) 
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Indeed, by the independence of 


E(n-iT„)2 = 

n n 

2=1 2 = 1 

Then, it follows from the bonndedness of moments in fl3.3p that 


pml. 


Ein-^T^y = 0{^^) + 0{ 


. {n — pn' + rnkY 


n 


n 


—)■ 0 , 


(3.15) 


which yields fl3.14p . as claimed. 


In view of (I3.14p and pn'/n —)■ 1, the proof of Theorem 13.21 now rednces 
to proving that 


ypn' ^ 

^ 2=1 


(3.16) 


For this pnrpose, let us hrst treat the limit behavior of the variance. Using 
the notations as in [ 23 ], we set 


rrik 


A . — WV^ -I- 9 \ ^ FV . V ■ 

z^n,i ■ n,i+mk ' ^ n,i+mk ^ n,i+mk-j ■ 

i=i 


(3.17) 


By the independence of {Un,i}i<i<p and the m^-dependence of j}i<j<„, 
straightforward computations show that 


V Un,iy = E(u„ 

ypn' pn' 




pn ^ 
i=l 


/ mfe ^ 

^ n'—27nk 

IE j ^ ^ — 


\h=l ) 

h=l 


2 1 ^ n'- 2 mk 

^(~r) - ^n,{i-l)n'+h, 

n' pn' 


(3.18) 


2=1 h=l 

where the last step is due to fl3.3p . 
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Set Ui^h '■= {i — + h + uik- In order to take the limit p, n' —)■ oo, we 

note from fl3.17p and the definition of that 


^ p n'-2mk 

^ X] X] ^n,{i-l)n'+h 
^ i=l h=l 


p n'-2mk 

-—Y y 

rm' 

^ i=l h=l 


ruk 


n,ni^h 


2 '^Yn,ni^h'^n,ni^ 

i=i 


h—J 


p n'—2mk 

(^r’^Vari 


pn 


i=l h=l 


n: 


ak 

i,h 


n rrik p n'-2mk 
^ j=i i=i h=i 


^k,ni h ^k,ni ^—j 


Yh ’ Yh - j)“^ 
It is not difficult to see that, for every 0 < j < rrik, as p, n' —)■ oo, 

1 


), (3.19) 


p n'-2mk 

/ '^i,h sak^ '^i,h ~ J \ak ^ 
pn' ^ ^ n n 

^ i=l h=l 

Moreover, we have that as n(f, h) —)■ oo. 


2ak + 1 


(3.20) 


f 3 \ ^ { ry ry \ 

7-^ ^ Cov{Zk,x, ^fc,i+ih 

Kh {^i,h-3) 


(3.21) 


where 1 < f + 1, are dehned as in fl3.6p . (See the Appendix for 

the proof.) 

Thus, plugging fl3.20p and fl3.2ip into fl3.19p . we obtain that as p, n! —)■ oo. 


Y p n'-2mk 

^ Y^ Y^ ^n,{i-l)n'+h —t Dk 
n i=i h=l 

with Dk defined as in (13.5p . which together with (I3.18p yields that 


(3.22) 




—>■ Dk. 


(3.23) 
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As regards the third moments of Un,i /in fl3.16p . since E|y„ 
and n' = o{p'3), we have for 1 < i < p, 




3 

in'—m]^ 

E 

1 — 
/—-Un,i 

yn' 

= (n')"^E 

j={i—l)n'-\-l 




{n' - mkf) 



=0{{n')h 

= o{p^). 


= 0 ( 1 ) 


(3.24) 


Consequently, combining fl3.23p and (13.241) . and applying Lyapunov’s cen¬ 
tral limit theorem (see [211 Appendix], see also [9l Theorem 7.1.2]), we obtain 
fl3.16p and complete the proof of the hrst part (i). 


{ii). In order to prove (13.71) . we just need to show that for all zi, • • • , Zr & 

M, the linear combination is normally distributed 

in the limit. 

To this end, as in the proof of (13.131) . the limit behavior is the same as 
that of 

r n nr 

n~^ ^ ^ Xk.^i = n~^ ^ ^ Zjn~°‘^^Xk^^i, (3.25) 

j=l i=l 2=1 j=l 

where Xk^^i is dehned as in (I3.12p with kj replacing k. 

Since ,*}*>! is m^-dependent for each 1 < j <r, Xk^^i}i>i 

is M-dependent with M = maxi<j<r . Thus, as in the preceding proof of 
the hrst assertion (i), we deduce that ^kj,i is nor¬ 

mally distributed in the limit, which consequently implies the Gaussian huc- 

tuation of the random vectors n~^{n~°'^^TrQ^, ■ ■ ■ ,n n>l. 


It remains to compute the covariances Let rriij = ma.x{mki,mkj}- 

Then, {Xk^^q, Xkj,q}q>i is m^-dependent. Set := n~‘^^*Xki^q, and 






h=l 
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Straightforward computations show that 


n 


=n 


— {a(ki+k 




kj,q) 


q=l q=l 




q=l 


q=l 


9=1 


n,q+mij Yk„ n,q+mij 


9=1 


mu n—m 


ij If' ''f'lj 


'kl ^ ^ ^ ^ {^^ki,n,q+mij—hy^kj,n,q+mij Y'^^ki,n,q+mijYkj ,n,q+rrnj—h^ ■ 

h=l q=l 

(3.26) 

As in the proof of fl3.22p . we deduce that for each 0 < h < rriij, as n —)■ cxo, 

n—mij 

^ q + rriij - h ^ g + rriij ^^ 1 

n ^ n n 

9=1 


oi(ki -\- kj^ -\- 1 


and as q ^ oo, 


Cov{ . ^ ^ Xk,,q+m,, Cov{Zk. 1 , Zk. l+h). 

{q + rriij — {q + rriijY^i ^ ’ ■” 


Hence, 


n 


'y ^ ^O^ki,n,q+mij-hYkj,n,q+mij) 


9=1 

Similarly, 


<y.{ki + /cj) + 1 




(3.27) 


n 


y ^ IE(^fei,n,q+mij ) h) t 


9=1 


a(/cj + /cj) + 1 


Cov(^Zk^^l+h) Zkj^l)- 


(3.28) 


Therefore, it follows from fl3.26p - fl3.28p that, as n ^ oo, 


^-{a(ki+kj, 


•^+»Cav(J2Xk„„Y,Xti.,) ^ A(i,j), 

q=l q=l 

where A{i,j) is dehned as in fl3.8p . thereby completing the proof. □ 
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3.3 Proof of Theorem 13.3 


Proof of Theorem \3.3[ The proof is similar to that of Theorem 13.21 we only 
need to prove that, with the additional scaling the limits fl3.9p and fl3.10p 
exist. Below, we show the existence of the limit fl3.9p . the proof of (13.1 Op 
follows analogously. 

Taking into account the dehnition of j in fl3.12p . we need to show the 
asymptotic estimate of 11^=0 ^n+j- 

Indeed, by Assumption {H.2), for each 0 < j < / — 1, 


=ab + + O(^) 


Here, with a slight abuse of notation, 0{n stands for the term of order 
after taking the expectation. Hence, for > 1, 





+ 





Hence, setting |ht; 


^ mj, we derive that for |??t/| > 1, 
i=i 


Qn+j'Pn+j ^m,+i 

i=o 


i=i i=i 


(3.29) 


Similarly, setting \n i\ = '^nj, we have that for \ni\ > 1, 

j=0 


tl(^r‘ - tl(d"‘+<1 


j=o 


j=0 


7-1.^ 


On 


]+ 0 {—) 


n 




n 


3=0 


Cn+i 1 




n 


(3.30) 
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Thus, fl3.29p and fl3.30p yield that for > 1, \lti\ > 1, 


i-i 


j=o i=o 


^ abd + '^^(adrnj^n+j-i + bdrnjr]n+j-i 
i=i ^ 


+ abrijCn+j ]n ^ + abrioCnn 


+ 0{n-^^). 


Moreover, it is easy to see that, for \r^i\ =0, 


(3.31) 


= + !•%) +0(4), 


(3.32) 


and for \lti\ = 0, 

i-i 

/ ^n-\- 

-t2a 


n , (^n+jbn+j 
^ n2 

3=0 


= a 


6) 2 ^[ab + a ^ 


m 


i=i 


in+j- 

"? 




i=i 


n 


(3.33) 


Plugging fl3.12p . fl3.3ip - fl3.33p into the right hand side of fl3.9p . we conse¬ 
quently deduce that the limit fl3.9l) exists, which actually can be computed 
explicitly in terms of the covariances of and For simplicity, we omit 
it here. The proof is complete. □ 


3.4 Applications 

Let us start with the independent identically distributed (i.i.d.) case, which 
is a direct consequence of Theorem 13.21 . 

Corollary 3.5 (The case a = O.j In the symmetric case, assume that {a*} 
is a sequence of i.i.d. random variables, {4} satisfy that di = /(aj_i, Oj) with 
f a continuous function, or {di} are i.i.d random variables and independent 
of {oi}. In the non-symmetric case, assume that {(aj_i, 4, 4)} is a sequence 
of i.i.d random vectors. Suppose also that in both cases the entries have all 
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moments finite. Then, the assertions in Theorem \3.2\ hold. In particular, the 
variance is given by 

Dk = Var{Xk^2) + 2 Cov{Xk^2, Xk^2+j), (3.34) 

i=i 

with mk defined as in Theorem \3.3\ and Xk^i as in fl3.12|] |^. 2 < i < mk + 2, 
and the covariances A{i,j) are given by 

rriij 

^ihj) = Cov{Xk^^2, Xkj,2) + '^^iCov{Xk^^2, Xkj,2+h) + C0v{Xk^^2+h, Xk.^2)), 

h=l 

where mij = max{mfc., 

Corollary 13.51 have several applications to the physical models and ran¬ 
dom matrix models mentioned in Section [TJ 

Anderson model m W' When Qi = bi = —1, and {di} are i.i.d 
random variables with all moments finite, this tridiagonal random matrix 
dH]) is referred to the Anderson model in the physical literature, restricted 
to the bound domain A„ = {1, • • • ,n}. Write Qn = Qn,o + Ki, where 
Vn = diag{di, ■ ■ ■ ,dn). Qn,o represents the finite difference operator, a dis¬ 
crete analogue of the one dimensional Schrodinger operator, and W repre¬ 
sents the operator of multiplication by the random field di, 1 < i < n. For 
this model, the conditions in Corollary 13.51 are verified. Hence, its traces of 
powers TrQ^, k > 1, are normally distributed in the limit. 

Hatano-Nelson model When Oi/hi > 0 and {(aj_i, dj, 5*)} is a 

sequence of i.i.d. random vectors with all moments finite, this tridiagonal 
random matrix is motivated by the non-Hermitian quantum mechanics of 
Hatano and Nelson (see [21] and references therein). In this case, it follows 
from Corollary 13.51 that the traces are approximated normally distributed. 

Random birth-death Markov kernel This tridiagonal random ma¬ 

trix arises in the random walks with random environment in chain graphs. 
Given the chain graph G = (V, E), V = {1, ■■■ , n}, E = {{i,i), \i — 3\ < !}• 

^The subindices of Xk^i here starts with 2, since, e.g. in the case di = f[ai-i,ai), i > 1, 
the distribution of (ii(= /(0,ai)) is different from that of di = f{ai-i,ai), i>2. 
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Consider the random birth-death Markov kernel with the state space V 
and the entries satisfying ai,bi G (0,1], dj G [0,1], and aj_i + di + bi = 1. 
One concrete example is the random conductance model, in which case 
bi = ttn-i = 1, bi = 1 - tti-i = Ui^i+i/{Ui^i+i + Ui^i-i), and {Ui^i+i} are 
i.i.d. positive random variables. Another example is that bi = an-i = 1, 
bi = 1 — aj_i = Vi, and {Vi} are i.i.d. random variables on [0,1]. We refer 
to [5] for more detail discussions and for the study of the limiting spectral 
distribution of such model in the ergodic environment. Here, in the i.i.d. 
environment, i.e. {{ai-i,di,bi)} is a sequence of i.i.d. random vectors with 
all moments hnite, we deduce from Corollary 13.51 that the traces are asymp¬ 
totically normally distributed. 

Random birth-death Q matrix ({23^ )■ This model is motivated from the 
inhnitesimal generator of continuous-time Markov process in the chain graph. 
In this case, Oj, bi > 0 and di = — (aj_i -|- R). In [23] we proved the existence 
of the limiting spectral distribution in the strictly stationary ergodic case 
(see [231 Theorem 1]). Here, in the case that {oj} and {bi} are two sequences 
of i.i.d random variables, {a*} is independent of {bi} (or Oi = bi, i > 1), and 
they have all moments finite. Corollary 13.51 implies the Gaussian fluctuations 
of the traces. 

Now, we formulate and prove Corollary 13.61 below for the case a > 0, 
which is mainly motivated by the symmetric tridiagonal random matrix stud¬ 
ied in [29] and includes, in particular, the well-known /3-Hermite ensemble. 

For simplicity, we will regard the I -|- 1-dimensional vector iti in in 

fl2.3p as the vector in the larger space by just adding zeros to the 

remaining [|] — I coordinates, and for each 0 < g < [|], denotes the g-th 
normal basis in namely, the g-th coordinate of ~^q is 1 and the others 

are zeros. 

Corollary 3.6 (The case a > t].) Consider the symmetric tridiagonal ran¬ 
dom matrix (II.ip . i.e. Oi = bi, i > 1. Assume that in {H.l) all a* and di 
are independent. Assume also {H.2) with 0 < e < a and sup„>i < oo 

for any k > 1. Then, the assertions in Theorem \3.3\ hold. In particular, the 
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covariances A{i,j) are given by 


Hhj) 


qkj+kj h h ( \ 

a{ki+kj)+l-2e \ki/2) \kj/2) ’ 
ah+l^j-2yar{Q , , / fc^-l W kj-l \ 

a{ki+kj)+l-2a V(A:i-l)/2/ \{kj-l)/2) ’ 

0 , 


if ki, kj even; 


if e = a, ki, kj odd ; 
otherwise. 


(3.35) 


Remark 3.7 This result coincides with Theorem 3 in /^P|/FI under the con¬ 
ditions considered here. However, the proof presented below is quite different 
from that in JWf, it is analytic and mainly based on Theorem \3.3[ 


Proof. First, by the condition on dn, we note that in Assumption {H.2) 
d = 0, (n = n~^^~^^dn, and the limit C = 0 when e < a. Below we shall 
discuss the case k is even or odd respectively. 

When k is even, fl2.2p indicates that i\ := is also even. Since 

d = 0, from fl3.3ip - fl3.33p we see that the main contribution comes from the 
case [iti\ = 0, i.e. there is no loop in the circuit. Hence, it follows from 
(Km and (Km that 

(3.36) 


which by the independence of {Cn} implies that 


'^ki,kj (1 + h, 1 ) 

=4a^^+^^-Var(0 




^ r^' it' 

"'ll 


'y ^ '^q^qidq+h,q'■ 

l<q<l 


(3.37) 


We shall compute explicitly the covariance A{i,j) in fl3.1ip . For conve¬ 
nience, we place the leftmost vertex of the circuit tt' of type {f 
at 1. Then, the circuit tt of the type {l,r^i,lti) contributing to fl3.37p is 
the one with the leftmost vertex 1 -|- h, and Ylii<q'<v '^q'^'q'dq+h,q' = 


Hn [2 
|(fc+0 


(3.5)], , ^Ij.(7i)+I/.( 72) denominator a(k 1) shall be modified by 


,5(Id7l)+Ii(72)) 


and a{k -|- ^) -|- 1 — 2a respectively. 
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l/4E,>o Ih+q{n)Ih+q{'^'), where Ij( 7 ) dentes the number of edges {i,i + 1 } 
in the circuit 7 , i.e. Ij( 7 ) = if the leftmost vertex of 7 is j (see 

also [29l p.l87]). Hence, the summation in + h,l) consists of the 

overlapped parts of all these two circuits vr and vr'. Similar argument holds 
for 1 + h), with the slight modihcation that the leftmost vertex of 

tt' is 1 — h (we here extend the positions of vertices to the negative in¬ 
tegers). Thus, using the notation T{k,l) in [221 (3.1)], we deduce that 
the summation in (13.lip is eaual to - )cr(’- ) Ei<oIi(7l)I*(72) = 

(see [29l p.213] for the equality), which together with the 

coefficients {a{ki + kj) -|- 1 — 2e)~^ and in (13.lip and (I3.37P 

respectively yields (I3.35P when ki, kj are even. 

When k is odd, ( 12 . 2 p implies that |7t/| shall be also odd. As in the 
preceding arguments, (I3.3ip - (l3.33p yield that the only case contributing to 
the main order is that \^i\ = 1 , namely, there is only one loop in the circuit. 
We use the g-th normal basis q of M[ 21 +^ to indicate that the vertex of this 
loop is n -|- g if the leftmost vertex of the circuit is n, 0 < g < [|]. Then, 
(I3.3ip yields that 

n-^^Xk,n= Y. + 0{n-^^). (3.38) 

o<g<[|] 


(Note that, this includes the easy case k = 1.) Hence, it follows that 

'^ki,kj (1 T h, 1) 


=a 


—2 


Var{C) Y Y 


X q ll X q' r 

'^l Oq+h,q'- 


(3.39) 


0<g<[|] 0<5'<[|] 


In particular, cr^. ^.(1 -|- h, 1) vanishes when 0 < e < a, since in that case 

C = o. 

The computation of A{i,j) in this case is easier. Similarly, place the left¬ 
most vertex of tt' of the type (/', at 1. (13.391) indicates that leftmost 

vertex of vr of the type (/, 7^^) is 1 -|- h, and the summation in (I3.39P con¬ 

sists of all these two circuits tt and tt' with exactly one loop at the same vertex. 
Thus, the summation in (13.lip is equal to \T{ki, kj)\ = {{k^Y)/^'^ 

(see [221 Remark 4] for this equality), which yields consequently (I3.35p . 
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Finally, as all and Cn are independent, it follows from fl3.36p and 
fl3.38p that cTfc.^fc.(l + h, 1) = 0 in the other cases, which implies immedi¬ 
ately A{i,j) = 0. The proof is complete. □ 


f3-Hermite ensemble (,T7, This model has been widely studied in 

the literature. In this model, ai = bi, ai is distributed as iXi/s 

is the X distribution with i(3 degrees of freedom), (3 > 0, di is normally 
distributed as N{0,2/{3), and {ai,di} are all independent. One signihcant 
fact is that the density of the eigenvalue distribution is given by Cnp ni<j 1-^*“ 
Xj\^ exp{—f3 wifh Cnfs a normalization. Hence, it generalizes (3 = 

1, 2,4 in the classical Gaussian ensembles (i.e. GOE, GUE, GSE respectively) 
to continuous exponents /3 > 0, which are connected to lattice gas theory, 
Selberg-type integrals. Jack polynomials and so on. Now, using the fact that 
Xr — -^(0,1/2), we deduce that in Assumption {H.2) a = e = 1/2, 

a = 1, r; = A^(0, l/(2/9)), d = 0 and ( = N(0,2/(3). Thus, Gorollary fl3.6p 
implies Gaussian fluctuations of the traces and the covariances are given by 


A(bj) 


I kjkj I ki \/ kj \ 

^ ki-\-kj \ki/2) \kj/2) ’ 
i_kikj_/ ki-1 W kj-l \ 
/3 ki+kj \{ki-l)/2/ \{kj-l)/2j ’ 

0 , 


if ki, kj even; 

if e = a, ki, kj odd ; 

otherwise. 


which coincides with [29l Gorollary 20 and [m Theorem 1.2]. 


4 Deviations 

This section is devoted to the large deviation and moderate deviation princi¬ 
ples for the traces. Due to technical reasons, we will consider the tridiagonal 
random matrix (II.ip with entries satisfying Assumption (iJ.3) in the case 
a = 0 below. 

(H.3) {i). In the non-symmetric case, di, 6i)} are i.i.d. random vectors. 

Set 1/ := P o (oi, ^2,62)”^ £ ^(M^). 

In the symmetric case (oj = 6*, Vi > 1), {ai} is a sequence of i.i.d. 
random variables, and {di} satisfies: 

^In [29l Corollary 2], the terms in the even and odd cases shall be modified by 

kl(T^ 1 feZcTj 1 

a{k+i)+i- 2 e and ^^k+i)+i- 2 c. respectively. 
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(a), {di} is a sequence of i.i.d. random variables and independent of 
{aj}. Set u' :=F o (di, ai)“^ G Or, 

(in), di = /(oj-i, tti), where / is a continuous function and oq = 0. Set 
z/" = Po (ai)-^ e ^(M). 

Moreover, in all cases {oj, di, bi} are bounded random variables. 

Remark 4.1 Assumption {H.3) already includes the models corresponding 
to the case a = in Suhsection \3.4\ such as the Anderson model, the Hatano- 
Nelson model, the random birth-death Markov kernel and the random birth- 
death Q matrix. 


4.1 Large deviations 

Let us first introduce some basic notations and definitions for large deviation 
principles. More details can be found in [13] and [20]. Then, we formulate 
and prove our result in Theorem 14.21 below. 

Consider a complete separable metric space ^denotes all the 
probability measures in A sequence of probability measures {pn} C 
is said to satisfy the large deviation principle (LDP for short) with 
speed Sn —)■ cxo and good rate function / : ^ —)■ [0,cxo], if the level sets 
{x & ^ I{x) < c} are compact for all c G [0, cxo) and if for all Borel set A 
of St:, 

1 1 

— inf I{x) <liminf — log/i„(A) < limsup — log/i„(A) < — in£/(x), 

x&A° n^oo Sn n^oo Sn xGA 


where A° and A denote the interior and closure of A respectively. In that 
case, we shall simply say that {fin} satishes the LDP{sn, I) on . We also 
say that a family of .^-valued random variables satishes the LDP{sn, I) if 
the family of their laws does. For any fi, fi' G the relative entropy of 

fi with respect to fS is dehned by 




j g log gd/i', if 9 = ^ exists; 
oo, otherwise. 


We also consider the product space SP^ with r > 2. A probability measure 
fi G is said to be shift invariant if for any Borel set A of 


fi{x G : {xi, ■ ■ ■ , Xr-i) G A) = fi(x G : (x2, • • • ,Xr) G A). 
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Moreover, for any fi G P ^ define the probability measure 

p®r P ^ by, for any Borel set A of 


{^i®rp) (A) = / n{dx) / I{ix,y)eA}p{dy), 


s:^ 


X 


where I[(^x,y)&A} is the indicator function of the set {{x,y) G A}. Dehne for 
p G ^(^), 


dr,p{p) 


H{p\pr-i ®r p), if P is shift invariant; 
oo, otherwise. 


(4.1) 


where Pr-i denotes the marginal of p on the hrst (r — 1) coordinates. 

Now, come to the tridiagonal random matrix fll.ip . For each /c > 1, set 
rfc = [|] + 1. In the non-symmetric case, take ^ = M^, and dehne the 
function F on by 


i-i i 

F(«i, • • •, Ylixi+jZi+jr^+^ n ^i+i’ 

t=o j=o 

where a* = {xi-i,yi, Zi) G R^, 1 < f < Tfc and xq = 0. Associating with F, 
dehne by 

■= '^^HAAp) ■ P e .^((R^)^"), X = piF)}, 

where u is dehned as in Assumption {H.3), and p{F) denotes the integration 
of F with respect to p. 

Similarly, in the symmetric case {ii) in {H.3), take ^ = R^, dehne the 
function F' on (R^)^* by 

Y1 ( 4 - 3 ) 

i=o i=o 

with f3i = (yi, Zi), 1 < i < r^, and the corresponding function I Ay is dehned 
by 

lAyA) := ini{Ir,,Ap) : p G ^((R^)^''), x = p{F')} 
where F is dehned as in Assumption {H.3). 
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In the symmetric case {Hi) in {H.3), take = M, and define the function 
F" on by 


F"(zo,-- - ,Zr^) = F'(7i,--- ,7^J, 

where 7 * = {f{zi_i,Zi),Zi), 1 < i < r^. The related function is 

defined analogously by 

Ci..”(^) := inf{t.+i,."(;i) : #1 e ^((R)’'>+'), a: = I,(F")} 
with u” defined as in {Hi) in {H.3). 

We are ready to state our large deviation result. Since when k = 1, 
TrQ^ = i-®- of random variables, of which the large 

deviation is well known by the Cramer theorem (see [131 Theorem 2.2.3]). 
Therefore, below we are concerned with the case k > 2. 

Theorem 4.2 Assume {H.3). Let k > 2 and set := [|] + 1. Then, 
satisfies the LDP{n, LDP{n, Ifi^^y) and LDP{n, Ir”+iy,) 

in the cases {i) — {Hi) in {H.3) respectively. 

Proof. Let us first consider the non-symmetric case {i) in {H.3). We 
note that {-TrQ^} and {-Yll=i^k,i\ is exponentially equivalent, i.e., for 
any 5 > 0 , 


1 1 

lim -logP(-|TrQ^ - > 5) = -cx). (4.4) 

n^oo n n 

i=l 

(See the Appendix for the proof.) Thus, by the exponential equivalence 
theorem! [131 Theorem 4.2.13]), the proof reduces to proving the large devi¬ 
ation of i YA=l^k,i- 

Now, set the random vectors Oj = (aj_i, d*, 6 j) G M^, i > 1. By the 
independence and the identical distribution of {oj} in Assumption {H.3), 
we regard {oi} as a Markov chain on with the transition probability 
7i{x,dy) := ^{dy), x,y E Note that {7i{x,dy)} obviously satisfy the 
uniform Assumption {U) in [131 P-275] (see also Hypothesis 1.1(a) in [20]L 
i.e., {cTj} is a uniform Markov chain on Thus, applying [131 Theorem 
6.5.12] (see also [20l Theorem 1.4]), we have that the multivariate empirical 
measure ya„ := ^ satisfies the LDP{n, with the good 

rate function defined as in fl4.ip . 
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On the other hand, by the definition of j in fl3.12p . we note that 

^k,i ' ' ' 1 —l)) (‘^•^) 

namely, j can be viewed as a continuous function on the product space 
Moreover, 

^ n 1 ^ 

^ ^ Xfc j ^ ^ F(^CXi, ■ ■ ■ , Oj-i-r’j. — 1) /^n(X), 

i=l i=l 

which implies that - 'Ym=i is an additive functional of the uniform Markov 
chain. 

Therefore, we can apply the contraction principle (see [131 Theorem 
4.2.1]) to obtain the large deviation of ^ as specified in the non- 

symmetric case. 

The symmetric case can be treated analogously. In fact, in the case 
(a) in {H.3), Xk^i = - and the random vectors (3i = 

{di, bi) G forms a uniform Markov chain on with the transition prob¬ 
ability 7i'{x,dy) = F{dy). Moreover, in the case (Hi) in {H.3), = 

F"{bi-i, ■ ■ ■ ,6j+rj,_i), and {bi} C M forms a uniform Markov chain in M 
with 7r"{x, dy) = u"{dy). Therefore, applying [T3l Theorem 6.5.12] and the 
contraction principle, we obtain the asserted large deviation results. The 
proof of Theorem 14.21 is complete. □ 


4.2 Moderate deviations 

Moderate deviation principles for dependent random variables are widely 
studied in the literature, see e.g. [3211II1E7] and references therein. Here, 
for the moderate deviations of the traces, we prefer to give an elementary 
proof based on the blocking arguments as those in the proof of Theorem 13.21 

Theorem 4.3 Assume {H.3). Set Xn = n~'^ with v G (0,1), and letTrQ^ = 
TrQ^ — ETrQ^. Then, for every k>l and > 0, where is defined as 

in (13.341) . {^J^TrQ^}n>i satisfies the LDP{X-^, ^). 


Proof of Theorem 14-31 Let X^^j = X^^j — EX^.^,. We first note that 
^TrQ^} and {\ is exponentially equivalent, i.e. for any 
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5 > 0 , 


limsup A^logP ( y — |TrQ^ - V'Xfc,*! > (5 I =-cx). (4.6) 

n ^ J 

(The proof is postponed to the Appendix.) Hence, it is equivalent to prove 
the moderate deviation for For this purpose, we will use the 

blocking arguments as in the proof of Theorem 13.21 

Let p = where £ G (|(1 — z/), 1 — z/). Set n' = [^] and r = n — pn'. 

Let Yn,i, Un,i, Zn^i and T„ be as in the proof of Theorem 13.21 but with a = 0. 
Then 



Moreover, denote by A„ the logarithmic moment generating function of 
i-e. An{t) = logEexp(fW^X](L^ Un,i), t G M. 

We shall prove below that, for every h > 0, 


lim sup A„ log P 

71^00 



— (X), 


(4.7) 


and 

lim A„A„(A”^f) = ^—Dk, (4.8) 

n—^oo Z 

where Dk is the variance defined as in fl3.34p . Then, by the exponential 
equivalence, fl4.7p implies that we only need to consider the moderate devi¬ 
ations of Un,i- Consequently, fl4.8p and the Gartner-Ellis theorem (see 
e.g. [I3]) yield the asserted moderate deviation principle for the traces. 


It remains to prove fl4.7p and fl4.8p . For the proof of 04.71) . since T„ = 
^n,i and Zn,i are independent, using Assumption (iL.3) and the Bern¬ 
stein inequality (see e.g. [21 p.21]), we have 

F{\\l^fn\ > 6) < 2e-^"A2{B^n+cS))^ 

V n 
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where Bl = ^E|TJ2 = 0{^) and c = sup ||Z„,,|U = 0{.[^). 

V l<i<p ^ 

Thus, 


\n\ogn\\l^Tn\ > 5) 
V n 


=Anlog2 


^(f + 


—oo, 


/\nn> 


which implies fl4.7p . as claimed. 


Coming to the proof of fl4.8p . we note that by the uniform boundedness 
in Assumption {H.3) and e > \{1 — v), 

Un, = 0{n'/,/^) = = o(l). 


Hence, it follows that 


A„A„(A„^t) ='^ An logEexp(^=74,i) 
^ VAnU 

p 

= ^ An log 

i=l 

Similarly to flS.lSp and fl3.19p . we have 


B ~ EU^ ■ 

^ + ‘^<774' 

[Ann)2 


^ EUL ^ 


XnU 


/ "Ifc ^ 

^ n'—‘IvTi]^ 

IE 1 ^ Yn,{i—l)n'+h 


\h=l ) 

h=l 


XnU 


^ ^ n'- 2 m.k 

^n.P-lK+h, 

h=l 


(4.9) 


where An,(i_i)n'+ft are dehned as in (I3.17p . Hence, it follows from Assumption 
(7/.3) and similar computations as in fl3.19p - fl3.23p that 




XnU 


(4.10) 
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Moreover, for the third moment in fl4.9p . by the choices of A„ and n', 



in'—mj^ 


(4.11) 


Consequently, plugging fId.lOp and fld.lip into fl4.9p . we obtain that 



which implies fl4.8p . thereby completing the proof. 


□ 


5 Discussions 

1. For more general finite diagonal random matrix (see e.g. 0. H). the 
Gaussian fluctuations and deviations of the traces can be treated in a similar 
way. Indeed, using the finite width of band, we can still reduce the asymptot¬ 


ical analysis of the traces to those of the corresponding m-dependent random 
variables. 


2. It seems difficult to prove deviation results for the general case a > 0. 
In some special cases, e.g. the CUE case (a = 1/2), deviation results are 
known for the empirical spectral distribution (hence for the traces by the 
contraction principle), we refer to [Tl [T^. 

3. In the derivation of the large deviation results in Subsection 14.11 we 
regard the trace as an additive functional of a uniform Markov chain. With 
this point of view, one can expect to achieve similar results for the traces 
when the entries of fll.ip form an appropriate Markov chain, e.g. positive 
Harris recurrent. We refer the interested reader to [10] and references therein. 
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6 Appendix 


Proof of fl3.13p . By Assumption (H.l) and {H.2), 

supE|n-“^Qi^'’^f < cx), 

i>l 

where = off - EqJ. Hence, 


\n 


, —Oik 


{TrQi 


2=1 


Cp^‘ Y1 ri-'^'‘QtP'h 

i=n-l+l 


< Y1 

i=n-l+l 






with 11-112 denoting the standard norm, which implies fl3.13p . 


( 6 . 1 ) 


□ 


Proof of fl3.2ip . It is equivalent to prove that, as n —)> oo, 

n~‘^°‘^Cov{Xk,n, Xk^n+j) -t Cou(^fc,i, Zk,i+j). (6.2) 

First consider the nonsymmetic case. Let a* = (aj_i, di, bi), i > 1 and Oq = 0. 
By the independence and weak convergence of in Assumption {H.l), it is 
not difficult to deduce that 

(o^n; ■ ■ ■ ) t (cTl, ■ ■ ■ ) Cli_|_[^])) 

where 5j are independent but with the same distribution as that of {a + r],d + 
C,b + f). Then, by fl4.5p and the continuous mapping theorem ([191 Theorem 
3.2.4]), 

4-F(5i, • • • ,5i_^[fc]) (6.3) 

with F the continuous function dehned as in fl4.2p . 

Similarly, 

^ Xk,n+j —t F{ai+j, • • • , (6.4) 
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On the other hand, by 03.31) and Holder’s inequality 


supE(n < oo, (6.5) 

n>l 

which implies the uniform integrabilities of n~‘^°‘'^Xk^nXk,n+j, and 

rr°'^Xk^n+j, n>l. 

Therefore, by 06.3p - 06.5p . we can apply the Skorohod representation the¬ 
orem and the uniform integrability to take the limit and obtain 06.2p for the 
non-symmetric case. (See e.g. [T^ Theorem 3.2.4] for similar arguments for 
the bounded continuous mapping.) 


The symmetric case can be proved analogously. In fact, with the uniform 
integrability 06.5p . we only need to check the weak convergence of 

In the case that all di are independent of ai{= hi), let fdi = (dj, ai), i > 1. 
In this case, X^^n = F'ifdn, ■ ■ ■ ,/9n+[^]) F' dehnes as in 04.31) . Then, it 
follows from similar arguments as above that 

n-“'’X,,„4F'(;gi,--- ,^,+[|]), (6.6) 

where A, 1 < f < [|] -l- 1, are independent but with the common distribution 
as that of (d a -|- r^). 

In the case that d* = f{ai-i,ai), then /di = (/(ai_i, a*), a*), Xk,n is 
now a continuous function of a„_i, • • • ,a^^^ky As (a„_i, • • • ,a^_^_^k^) -A 
(oi, • • • by the continuous mapping theorem, we can obtain the 

weak convergence of n~^^Xk^n- The proof is consequently complete. □ 


Proof of 04.4p . Set N{k) = |Tfe|, the number of sets in 4/^ which is defined 
as in 02.31) . N{k) is finite and depends only on k. Note that 


P 



TrQi 

i=l 



< E p 


QpLi,iii 


n 

E 


i=n—l-\-l 



Then, letting C denote the maximum of 


a 


rii 


l,ni 


over the finite sets in 4/^ and 
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setting 6k = 5/([|]CiV(/c)), we deduce that 


P 



TrQi-Y^Xt, 

i=l 



< 


E E *• 

i=n-l+l 


Qi‘ 



Since by fl2.ip and Assumption {H.3), {IQh*’""*|} i>i is uniformly bounded, 

which implies that > nSk) = 0 for n large enough. Hence, by 

Lemma 1.2.15 in [13], 


-logp(-|Trg^->5 
n \n 


2=1 


< ,, max i logP (\Qi’%\ > n5k 


= —cxo. 


yielding fl4.4p as claimed. 


(6.7) 

□ 


Proof of (14.6p . Similarly to the proof of (16.7p . we derive that 


An log P 



n 


|TrQ^-5^Xfc,|><5 

2=1 









with 6k dehned as in the proof of (16.71) . which yields (14.61) . due to the fact 
that is uniformly bounded and n/Xn oo. 
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